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Abstract 

We introduce a deterministic model for scale- free networks, whose degree 

distribution follows a power-law with an exponent 7. At each time step, each 

vertex generates its offsprings, whose number is proportional to the degree 

of that vertex with proportionality constant m — 1 (m > 1). We consider 

the two cases: first, each offspring is connected to its parent vertex only, 

forming a tree structure, and secondly, it is connected to both its parent and 

grandparent vertices, forming a loop structure. We find that both models 

exhibit power-law behaviors in their degree distributions with the exponent 

7 = 1 + ln(2m — l)/lnm. Thus, by tuning m, the degree exponent can 

be adjusted in the range, 2 < 7 < 3. We also solve analytically a mean 

shortest-path distance d between two vertices for the tree structure, showing 

the small-world behavior, that is, d ~ lniV/ln/c, where N is system size, 

and k is the mean degree. Finally, we consider the case that the number of 

offsprings is the same for all vertices, and find that the degree distribution 

exhibits an exponential-decay behavior. 
PACS numbers: 89.70,+c, 89.75.-k., 05.10.-a 
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I. INTRODUCTION 



Recently, complex systems have received considerable attention as an interdisciplinary 
subject [0,0. Complex systems consist of many constituents such as individuals, substrates, 
and companies in social, biological, and economic systems, respectively, showing cooperative 
phenomena between constituents through diverse interactions and adaptations to the pat- 
tern they create MM . Recently, there have been a lot of efforts to understand such complex 
systems in terms of networks, composed of vertices and edges, where vertices (edges) repre- 
sent constituents (their interactions). This approach was initiated by Erdos and Renyi (ER) 
||. In the ER model, the number of vertices is fixed, while edges connecting one vertex 
to another occur randomly with a certain probability. However, the ER model is too ran- 
dom to describe real complex systems. Recently, Barabasi and Albert (BA) |||7|] introduced 
an evolving network where the number of vertices N increases linearly with time rather 
than fixed, and a newly born vertex is connected to already existing vertices, following the 
so-called preferential attachment (PA) rule; When the number of edges k incident upon a 
vertex is called the degree of the vertex, the PA rule means that the probability Ilj for the 
new vertex to connect to an already existing vertex i is proportional to the degree hi of the 
selected vertex, that is, 

The main difference between the ER and BA models appears in the degree distribution. 
For the ER network, the degree distribution follows the Poisson distribution, while for the 
BA network, it follows a power-law, P(k) ~ k~"< with 7 = 3. The network whose degree 
distribution follows a power-law is called the scale-free (SF) network 0. SF networks are 
abundant in real- world such as the world-wide web @-(llJ, the Internet [ 12 -pL5 ] , the citation 



network fl6|| , the author collaboration network of scientific papers [fuH , and the metabolic 
networks in biological organisms ||18|. 
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While a lot of models have been introduced to describe SF networks in real world, most 
of them are stochastic models. However, a couple of models recently introduced by Barabasi, 
Ravasz and Vicsek (BRV) |19[, and Dorogovtsev and Mendes (DM) |2| are deterministic. In 
general, the deterministic model is useful in investigating analytically not only topological 
features of SF networks in detail, but also dynamical problems on the networks. Both the 
BRV and the DM models are meaningful since they are not only the first attempts for deter- 
ministic SF networks, but also the ones constructed in a hierarchical way, so that analytic 
treatments can be made easily using recursive relations derived from two structures in suc- 
cessive generations. In the BRV model, however, the mean shortest-path distance between 
two vertices, called the diameter, is independent of system size. Thus, the model may be 
relevant to some specific systems such as the metabolic network [fDjj] , where the diameter is 
independent of system size. In this paper, we introduce another type of the deterministic 
model for the SF network, which is also constructed in a hierarchical way. Our model is 
based on almost the same idea as that of the DM model. While the DM model starts from 
a triangle, our model does from a tree structure. This difference makes one easily modify 
the model into more general cases such as loopless or loop structures, and the ones with 
a various number of branches. Moreover, the simplicity of our model enables us to obtain 
the analytic solution for the degree distribution and the diameter. In particular, our model 
includes a control parameter, so that by tuning the parameter, we can obtain SF networks 
with a variety of degree exponents in the range, 2 < 7 < 3. Therefore our model should be 
useful to represent various SF networks in real world. 

This paper is organized as follows. In section II, we will introduce deterministic models 
specifically for tree and loop structures, respectively. In section III, analytic treatment will 
be performed for the deterministic models introduced in section II. The final section will be 
devoted to conclusions and discussions. 
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II. DETERMINISTIC MODEL 



It is known that the number of vertices in most of SF networks in real world increases 
exponentially in time. Thus, our deterministic model is constructed in an evolving way, 
where each already existing vertex produces its offsprings, and the connections are made 
between old and new vertices. Thus, vertices are generated in a hierarchical order, so that 
the number of vertices increases geometrically in time. On the other hand, it is known 
PD| , |2TH that the PA probability IT^, Eq.(|I|) is generalized for real networks as 

n , = kj + fi ^ 

y^. j kj -\- fi 

where fi accounts for some randomness in connecting edges. To take into account of this 
modified PA behavior, we introduce two rules, called the addition and the multiplication 
rule, in the deterministic model, depending on how new vertices are generated from each 
old vertex. The details on both rules will be described below. 



A. Tree structure 

The network forms a tree structure when new vertices generated from an old vertex are 
connected to their parent only. 



1. The addition rule 

In the case of the addition rule, at each time step, a constant number of new vertices, say 
i new vertices, are generated from each already existing vertex, and they are connected to 
their parent only. Then the degree k iA at vertex i, where the subscript a means the addition 
rule, evolves as 

k ha (t + 1) = k ija (t) + £, (3) 

so that 
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h, a (t) = i + e(t-ti), (4) 

for t > ti, where ti means the time when the vertex % was born. Then the total number of 
vertices newly born at time t becomes C a (t) — £(1 + €f for t > 1, and the total number of 
vertices N a (t) present at time t is 

N a {t) = Y l c a {j) = {i+er\ (5) 

j=0 

where N a (0) = 1 + £ is chosen. The definition of this model is illustrated schematically in 
Fie.l. 



2. The multiplication rule 

In the case of the multiplication rule, the number of offsprings generated from each old 
vertex is not the same, but it depends on the degree of each vertex. Let ki jTn (t) be the 
degree of vertex i at time t, where the subscript m means the multiplication rule. Then the 
number of offsprings generated at time t + 1 from the vertex % is proportional to its degree 
at the previous time, i.e., (m — l)k itm (t), where m — 1 is a proportionality constant. Thus 
the degree of vertex % increases by a factor m for each time step, that is, 

h,m{t) = mki tm (t - 1). (6) 

Thus, the degree of vertex i at time t is 

k i , m = m t - t \ (7) 

for t > tj. The total number of vertices newly born at time t, C m (t) can be obtained to be 

C m (t) = 2(m - l)m £ (* ~ l \m*\m - l)*" 1 ^, 
p =o V P J 

= 2m(m-l)(2m-l) w , (8) 

for t > 1. The total number of vertices N m (t) at time t is given by 

t 

N m (t) = ]T C m (t) = 1 + m(2m - 1)'. (9) 

3=0 
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The definition of this model is illustrated schematically in Fig. 2. 

One may write the rate equation for the degree in this multiplicative process with con- 
tinuous time as 

dh 



(m-l)ki tm . (10) 



at 

It would be interesting to compare this rate equation with the one for the preferential 
attachment (PA), in which the degree ki of vertex % evolves as 

— = C ki (11) 

where £ m (t) means the total number of edges newly introduced at time t. Since the total 
number of the degree at time t is given by 

Y,kj = 2m(2m- 1)*~\ (12) 

3 

and C m is given by Eq.(j|), Eq.flllD is reduced to Eq. ([[(]), indicating that the rate equation 
for the degree in the multiplicative process is equivalent to the one for the preferential 
attachment. 



B. Loop structure 

The loop structure can be formed in networks, when a newly born vertex is connected 
to more than one existing vertices. For the loop structure, each already existing vertex 
generates the same number of offsprings as those for the tree structure. However, a newly 
born vertex is connected to two distinct old vertices: one is its parent, and the other is its 
grandparent. When the parent is one of vertices on m branches (the centered one) born at 
t — 0, the centered one (one of vertices on m branches in a symmetrical way) is regarded as 
a grandparent. This rule is valid for both cases of the addition rule and the multiplication 
rule. The details on the connection rule is illustrated in Fig. 3. 
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III. ANALYTIC SOLUTION 



A. The degree distribution for the tree structure 

Since the degree of a vertex has been obtained explicitly as in Eqs.(|J) and (|7|) and they 
are ordered with time, we can obtain the degree distribution using the relation, 

Pit) = qi-^(*.w>*)], (13) 

which is valid for both cases of the addition and the multiplication rules. The detail of 
analytic treatments for the degree distributions for both cases are given as follows. 

1. The addition rule 

Using the fact, P a ,t(h,a > k) — P a ,t{U < t — (k — l)/£), where the subscript t means the 
tree structure, we obtain that 

P t-(k-i)/e 

Pa,t(ha(t) > k) = E (1+4* 

= (i + £)- {k - 1)/e - (i + ey( t+1 \ (14) 

Applying Eq.(|13|) to Eq.([L4|), we obtain the degree distribution to be 

P a , t (k)<x(l + £)-^ e \ (15) 
So, the degree distribution P a t (k) in the addition rule decays exponentially with k. 

2. The multiplication rule 

Since the degree ki has been obtained explicitly as a function of time in Eq.(^), P m ,t{ki > 
k) is written as P m ,t(ki > k) — P m ,t(ti < t), where r — t — ln/c/lnm. Since P m ,t(ti < t) 
means the density of the vertices born earlier than r, 



7 



P m Ah>k) = ± £, " { '< ] 



=0 N m(t) 



I, 2(m- l)(2m- l)*'" 1 



ti=0 



+ m(2m — 1) 



K fc -ln(2m-l)/lnm_ ( lg ) 



Thus the degree distribution is obtained to be 

d[l - P m ,t(Ht) > 



Pm,t(k) 



dk 

oc k~^ m) , (17) 



where 

j(m) = 1 + ln(2m- l)/lnm. (18) 

In the limit of m — > 1, we get 7(1) = 3, while as m goes to infinity, we get 7(00) = 2. Thus 
by tuning the parameter m, we can get a variety of SF networks with different exponents in 
the range, 2 < 7 < 3. 



B. The degree distribution for the loop structure 

1. The addition rule 

Let rij j0 (t) be the degree of vertex i at time t for the loop structure, where a means the 
addition rule. Each old vertex receives edges from its I children and £ 2 grandchildren as 
they are born. So, Eq.(|3|) is modified as 

n i>a {t + l) =n i>a {t) + {£ + £ 2 ). (19) 

Thus, the degree distribution shows an exponential-decay behavior, 

P a j(n)oc(l + l + £ 2 y n/(e+e2 \ (20) 
where the subscript I means the loop structure. 
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2. The multiplication rule 

Let rii, m (t) be the degree of vertex i at time t in the multiplication rule for the loop 
structure. The degree of vertex i can be obtained, 

n^ m {t) = n i)Tn (t - 1) + (m - l)k i>m {t - 1) 

+ [m-lfh^it-2), (21) 

where the second term on the right hand side of the above equation results from the children 
of the vertex i, and the third term from the grandchildren of the vertex i. Thus, the degree 
at the vertex % becomes 

n itm (t) = 2m t ~ u - m M,_1 - m » ( m ~ )m t ~ u . (22) 

Since the degree n i m (t) depends on time t similarly to Eq.([7|), we can apply Eq.(|TJ) even to 
the loop case, except that r is replaced by r = t + ln(2m — l)/lnm — 1 — Inn/ In m. This 
replacement, however, does not affect the degree exponent at all. Thus, even for the loop 
structure, the degree exponent is reduced to the same value, 7 = 1 + ln(2m — l)/lnm in 
Eq.flD. 



C. The diameter for the tree structure 

The diameter d(t) is defined as a geodesic distance between two distinct vertices along 
the shortest path averaged over all pairs of vertices at time t, that is, 

d(t) = , w 1 r Y di i (t) , (23) 

where dij(t) is the shortest-path distance between vertex i to j. For simplicity, let T>(t) 
denote the sum of the shortest-path distances between two vertices over all pairs, that is, 

V{t) = Y,diM- (24) 

It is not easy to obtain a closed formula for T>(t) for both the tree and the loop structure, 
however, we list T>(t) for the tree structure in a few early times in Appendix. We trace the 
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formula for the tree structure in two limiting cases, m — > and m — > oo, as follows. 



Let us first consider the case of m — > 1. For this case, we denote m = 1 + e and e<l. 
The total number of nodes iV(t) at time t is given by 

N(t) = 1 + (1 + 2e)*(l + e) 

^2 + (2t + l)e + C(e 2 ). (25) 

Moreover, the sum of all shortest-path distances T>(t) becomes 

V(t) ^2 + 4(2t+l)e + C(e 2 ). (26) 

Using the relation Eqs.(|23|), we can obtain the average distance to be 

+ t (27) 



7 + 61og(A^-l) 3 
Therefore, the diameter converges to 4/3 in the limit of iV — » oo 



Next, we consider the case ofm-> oo. In this case, the term in the highest order of m 
could be dominant, so that we trace the coefficient of the term in the highest order of m as 
a function of time. 

P(0) = 2m 2 + lower order terms 

V(l) = [(2 + 3) + (3 + 4)]m 4 + lower order terms 

P(2) = [(2 + 2- 3 + 4) + 2- (3 + 2- 4 + 5) + (4 + 2- 5 + 6)]m 6 

-(-lower order terms 
p(3) = [(2 + 3 ■ 3 + 3 ■ 4 + 5) + 3(3 + 3 ■ 4 + 3 • 5 + 6) 

3(4 + 3- 5 + 3- 6 + 7) + (5 + 3- 6 + 3- 7 + 8)]m 8 

+lower order terms 
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^(0 = EQE(fc)(*+p+2K (t+1) 

+lower order terms. (28) 

Therefore the term in the highest order of m for V(t) is obtained explicitly to be 

V(t) « (t + 2)2 2t m 2 < t+1 >, (29) 

where the coefficient (£ + 2)2 2 ' means the number of pathways having the distance 2(t + 1), 
which is the farthest one at time t in the system. On the other hand, 

N(t)(N(t)-l)^2 2t m^ t+1 \ (30) 

Therefore, the diameter d(t) at time t becomes simply 



d(t) ^t + 2 

„ log(iV - 1) 



+ 2. (31) 



log(2m - 1) 

Thus, for large N, the above equation is rewritten simply as 

d(N) ~ In TV/ In A; (32) 
with the mean degree k ~ 2m, which confirms the small-world behavior. 

IV. CONCLUSIONS AND DISCUSSIONS 

We have introduced a deterministic model for the scale-free network, which is constructed 
in a hierarchical way. At each time step, each already existing vertex produces its offsprings, 
whose number is proportional to the degree of the vertex. Depending on whether each 
offspring is connected to only one or more than one old vertices, the network forms either a 
tree structure or a loop structure, respectively. We have obtained the analytic solution for 
the degree distribution and the diameter explicitly for the deterministic model. By tuning a 
control parameter in the model, we can adjust the degree exponent in the range, 2 < 7 < 3. 
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Thus this model can represent a variety of SF networks in real world. Moreover, we obtained 
the diameter of the deterministic model analytically to be d ~ In N/ Ink, where N is the 
system size and k is the mean degree. Since the network is generated in a hierarchical way it 
is expected that a variety of physical problems can be solved through this deterministic model 
by constructing recursive relations derived from two structures in successive generations. 
On the other hand, the deterministic model has a shortcoming that it does not include any 
long-ranged edge, connecting two vertices belonging to different branches separated at t — 0. 
Thus, while this model can be easily generalized by controlling the number of branches m, it 
is extremely vulnerable, and can be broken into pieces by a simple deletion of the centered 
vertex. Despite this shortcoming, we think that our deterministic model could offer a guide 
toward generating more realistic deterministic model for SF networks. 
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VI. APPENDIX 

The closed formula for the sum of the shortest-path distance between two vertices, V(t) 
are shown for < t < 2. 

V(t = 0) = J\f m + A/" ,o[l + 2(m - 1)], 

V{t = 1) = jV [m 2 + 2m(m - 1)] + jV , [m + 2m(m - 1) + 2(m - 1) + 3(m - l)(m - 1)] 

+ M,o[l + 2m 2 - 2(m + 1) + 2m + 3m(m - 1)] 

+ A/i,o,o[l + 2(m - 1) + 3(m 2 - 1) + 4(m - l) 2 ], 
V(t = 2) = 7V [m 3 + 2(m 2 - m){m - 1) + 2m(m - 1) + 2m(m 2 - m) + 3m(m - l)(m - 1)] 

+ N ,o[lm 2 + 2m 3 - 2 + 2(m - l)(m - 1) + 3{(m - l) 2 + (m 2 - m){m - 1) 
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+ (to - l)(m 2 - m)} + 4(m - l) 3 ] 

+ A/i,o[to + 2(m 3 — 1) + 3{(m(m — 1) — 1)(to — 1) + m{m 2 — m) + to(to — 1)} 
+ 4m(m - l) 2 ] 

+ M,o,o["i + 2(m 2 - 1) + 3{(m 3 - 1) + (m - l)(m - 2)} 

+4{(m 2 - to)(to - 1) + (m - l) 2 + (to - l)(m 2 - to)} + 5(m - l) 3 ] 
+ AT 2) o[l + 2(m 3 - 1) + 3{m 2 (m - 1) + m(m 2 - m)} + 4m(m - l) 2 ] 
+ A/" 2 ,i, [l + 2(m - 1) + 3(m 3 - 1) + 4{(m(m - 1) - l)(m - 1) + m(m 2 - 1)} 
+ 5m(m — l) 2 ] 

+ M,o,o[l + 2(m 2 - 1) + 3{(m 3 - 1) + (to - l)(m - 1)} 

+4{(m 2 - 1)(to - 1) + (to - l)(m 2 - to)} + 5(to - l) 3 ] 
+ AT 2>1A0 [1 + 2(to - 1) + 3(m 2 - 1) + 4{(m 3 - 1) + (to - 1)(to - 1 - 1)} 

+5{(m 2 - 1)(to - 1) + (to - l)(m 2 - 1)} + 6(to - l) 3 ], 

where 

AT = 1, 

■A/o,o = m, 

A/i,o = ^ 2 — m, 

M,o,o = m ( m - 1), 

■A/2,0 = m 3 — m 2 , 

A/2,1,0 = ("^ 2 - m) (to - 1), 

A/2,0,0 = m(m 2 - to), 

and 

A/2,1,0,0 = m ( m - - !)■ 

A/i,j means the number of the vertices denoted by A iy j in Fig.2, where the first index % stands 
for its birth time and the rest indices j do its parent vertex. 
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FIGURES 





t=0 



t=l 



FIG. 1. Tree structures in the addition rule with I = 3 at t = and t = 1. 
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t=0 t=l 

FIG. 2. Tree structures in the multiplication rule with m = 3 at t = and i = 1. Ao stands 
for the vertex at center, Ao,o, a neighbor of Ao born at t = 0, and A^o (Ai^^), an offspring of Ao 
(Aq,o) born at t = 1. 
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FIG. 3. Loop structure in the multiplication rule with m = 3 at t = 1. 
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